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Rose bush architecture, among other factors, such as plant health, determines plant visual 
quality. The commercial product is the individual plant and interplant variability may be 
high within a crop. Thus, both mean plant architecture and interplant variability should be 
studied. Expansion is an important feature of architecture, but it has been little studied at 
the level of individual organs in rose bushes. We investigated the expansion kinetics of 
primary shoot organs, to develop a model reproducing the organ expansion of real crops 
from non-destructive input variables. We took interplant variability in expansion kinetics 
and the model's ability to simulate this variability into account. Changes in leaflet and 
internode dimensions over thermal time were recorded for primary shoot expansion, on 83 
plants from three crops grown in different climatic conditions and densities. An empirical 
model was developed, to reproduce organ expansion kinetics for individual plants of a 
real crop of rose bush primary shoots. Leaflet or internode length was simulated as a 
logistic function of thermal time. The model was evaluated by cross-validation. We found 
that differences in leaflet or internode expansion kinetics between phytomer positions 
and between plants at a given phytomer position were due mostly to large differences in 
time of organ expansion and expansion rate, rather than differences in expansion duration. 
Thus, in the model, the parameters linked to expansion duration were predicted by values 
common to all plants, whereas variability in final size and organ expansion time was 
captured by input data. The model accurately simulated leaflet and internode expansion 
for individual plants (RMSEP = 7.3 and 10.2% of final length, respectively). Thus, this study 
defines the measurements required to simulate expansion and provides the first model 
simulating organ expansion in rosebush to capture interplant variability. 

Keywords: Rosa hybrida L, individual plant, phytomer, model, elongation kinetics, leaflet size, internode length, 
growth 



INTRODUCTION 

Plant architecture constitutes the interface by which the plant 
gathers resources and perceives signals from its environment, 
which in turn modify plant architecture. In ornamental crops, 
such as rose bush, plant architecture is important in its own right, 
because it conditions plant visual quality, largely accounting for 
consumer choice (Boumaza et al., 2010). The management of 
plant architecture, by manipulation of the environment, in par- 
ticular, is therefore an important issue for rose bush growers. 
The development of plant architecture should be studied at two 
levels: that of the crop, the level at which management oper- 
ates and interactions occur between neighbors, and the individual 
plant level, as it is individual plants that are sold. A knowledge 
of mean values for a crop is therefore not sufficient; the variabil- 
ity between plants within a crop should also be characterized. To 
investigate the relationships between rose bush architecture and 



its environment, a functional-structural plant model (FSPM) of 
rose bush would be a powerful tool. Firstly, FSPMs account for 
plant architecture at the individual plant level; secondly, when 
coupled with phylloclimate models, FSPMs enable to estimate the 
physical environment actually perceived by aerial organs, which is 
heterogeneous within and between plants (Chelle, 2005). 

Plant architecture firstly depends on bud fate, which deter- 
mines the number and location of shoots; furthermore on the 
initiation, expansion and (re)orientation of internodes and leaves, 
and on floral transition, which determines the time of flowering 
and the number of flowers. All these traits have a potential impact 
on the visual quality of rose bush. The architectural traits that we 
wished to model in this study relate to the expansion of the main 
organs of the aerial vegetative apparatus: the stems and leaves. In 
addition to their physiological functions, the stems determine the 
shape of the plant (top-sided shape, symmetry). Together with the 
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leaves, they determine the compactness of the plant. Both these 
traits are thus important elements of visual quality (Boumaza 
et al., 2009). The expansion of an individual organ is charac- 
terized by the time of organ expansion, the expansion rate and 
expansion duration. Within a shoot, these traits follow gradients 
according to organ position which need to be characterized. 

Expansion features vary between species, and specific stud- 
ies are therefore required. The kinetics of leaf and internode 
expansion at the phytomer level have never been described for 
rose. Several models describing architectural variables have been 
developed for cut-flower roses. A first group of models pre- 
dicts variables at the level of the whole shoot, such as total stem 
length and basal stem diameter (Hopper et al., 1994; Costa and 
Heuvelink, 2003; Oki et al., 2006) or morphological quality classes 
(Morisot, 1996). Structural descriptions have been refined further 
in more mechanistic models (Lieth and Pasian, 1991; Dayan et al., 
2004), mostly with a view to predicting the harvest of rose flowers 
more accurately. Considerable effort has thus gone into model- 
ing photosynthesis and assimilate partitioning, but the kinetics 
of expansion of individual organs and its variation with respect 
to organ position have not been described. The FSPM of Buck- 
Sorlin et al. (2011) for cut-flower roses included the expansion 
kinetics of individual organs, but the paper focuses on simulating 
the local light climate and photosynthesis and does not present 
results for organ expansion. In rose bush, Demotes-Mainard et al. 
(2009) have described the coordination of leaflet and internode 
expansion kinetics as a function of position along the primary 
shoot. However, these preliminary results were obtained for only 
one set of growth conditions and relate exclusively to the behav- 
ior of an average 1 1 -phytomer plant. With a view to developing a 
FSPM for rose bush one step consists in modeling organ expan- 
sion. Therefore, additional knowledge of the kinetics of individual 
organs must be obtained. 

As in many species, organ size in rose is influenced by envi- 
ronmental factors, such as water (Demotes-Mainard et al., 2013), 
or nitrogen (Ashok and Rengasamy, 2000; Huche-Thelier et al., 
2011) availability, light quality (Rajapakse and Kelly, 1994; Maas 
and Bakx, 1995) and intensity (Hopper and Hammer, 1991; 
Bredmose, 1993; Maas and Bakx, 1995), mechanical stimula- 
tion (Morel et al, 2012), and genotype (Morel et al, 2009). 
Current knowledge of the effects of environmental factors is not 
sufficient to predict internode or leaf size in a range of environ- 
ments. However, one way to reconstruct the plant architecture of 
experimental crops accurately for the investigation of plant func- 
tioning involves using a model that simulates architecture with 
data obtained from experimental crops. These experimental data 
must capture the variability induced by the environment and not 
predicted by the model. This approach has been applied success- 
fully in several studies [for example in Baccar et al. (2011), to 
study the effect of wheat architecture on Septoria tritici epidemics; 
Kahlen and Stiitzel (2011), to study the photo-modulation of 
cucumber internode elongation] . The input data used to recon- 
struct architecture are generally the means obtained for plant 
samples. We suggest that the use of data from individual plants 
should make it possible to explain and reproduce not only the 
mean plant characteristics, but also the interplant variability. 
Using individual data would not increase the time required for 



data acquisition, because architectural data are necessarily mea- 
sured at the individual plant and organ level. For the correct 
reproduction of individual plants, the input data should capture 
interplant variability. This requires good knowledge of variations 
of the expansion kinetics of individual organs between plants. 

The aim of this work was (i) to investigate the expansion kinet- 
ics of rose bush primary shoot internodes and leaves and its 
variation between phytomer positions and between plants and 
(ii) to propose an empirical model reproducing organ expansion 
kinetics that accounts for interplant variability and makes use of 
non-destructive, easy-to-measure input variables. 

MATERIALS AND METHODS 

PLANT MATERIAL AND GROWING CONDITIONS 

Rosa hybrida "Radrazz" rose bushes were grown in Angers, 
France, in three experiments. Experiment 1 took place in spring 
2007 (2007-Sp), experiment 2 in summer 2007 (2007-Su), and 
experiment 3 in spring 2010, with a low (2010-LD) or high 
(2010-HD) plant density. Single-node cuttings bearing a five- or 
seven-leaflet leaf were harvested from the medial part of mother 
plant stems and grown for 4-5 weeks in humid conditions until 
rooting was achieved. Well rooted cuttings were planted in indi- 
vidual pots containing a mixture of neutral peat, coconut fibers 
and perlite. Plants were transferred to a greenhouse before bud 
break and grown, with a border row, at a density of 23 plants 
mT 2 in experiments 1 and 2, and at a density of 21 or 100 
plants mT 2 in experiment 3. The bud from the cutting pro- 
duced the primary axis. After floral bud and last leaf appearance, 
a variable number of lateral buds burst and developed in sec- 
ondary shoots. These shoots were let to grow but not measured. 
Plants were subirrigated, with tensiometer monitoring to ensure 
an absence of water stress. Mineral nutrition was provided by fer- 
tigation (5.0 mM KNO3, 2.0 mM Ca(N0 3 ) 2 , 2.0 mM NH4NO3, 
2.0 mM KH2PO4, 2.0 mM MgS0 4 , 0.25 mM NaOH; trace ele- 
ments (Kanieltra 6-Fe,0.1 ml.l _1 ,Hydro Azote, Nanterre, France); 
pH 5.6; EC 1.77mS.cm _1 ). Air temperature was measured above 
the canopy, in a ventilated shelter, with a platinum sensor. Leaf 
temperature was measured in experiment 3 on the abaxial surface 
of leaves, at about the height of the apex, with copper-constantan 
thermocouples. Photosynthetically active radiation (PAR) was 
continuously measured above each canopy, with a line quantum 
sensor (LI-191 LI-COR, Lincoln, Neb. USA). It was also measured 
with line quantum sensors at both the base and the top of the 
canopy, at various locations, at solar noon, under a cloudy sky, 
on three dates during primary shoot expansion in experiment 
3, for the calculation of intercepted PAR. The climatic condi- 
tions prevailing during the period of primary shoot expansion are 
presented in Table 1. 

PLANT MEASUREMENTS 

Destructive measurements of leaf dimensions 

Groups of 62, 17, and 18 plants were selected at random from 
the 2007-Sp, 2010-LD, and 2010-HD crops, respectively, and 
destructively sampled at various times (5-11 times, according 
to the crop), beginning when the basal leaves were expand- 
ing and ending when all the leaves were fully expanded. We 
determined the number of leaflets per leaf for all leaves (197, 201, 
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Table 1 | Prevailing climatic conditions and PAR interception during shoot expansion, at high and low plant density, in each experiment. 



Experiment 2007 spring 2007 summer 2010 spring 



High plant density Low plant density 



Incident PAR (mol.m^.day- 1 ) 8.0 ±1.6 11.7 ±1.1 7.7 ±1.1 



Intercepted PAR (%) 



Air temperature (daily mean, °C) 20.2 ±1.9 21 .7 ±1.3 

Leaf temperature (daily mean, °C) 

Humidity (%) 57.1 ±6.5 75.6 ±4.5 

Each value is the mean ± standard deviation. 

199 leaves for 2007-Sp, 2010-LD, and 2010-HD, respectively), 
the length and width of all leaflets and leaf length (measured 
with a ruler in 2007-Sp and by image analysis (ImageJ 1.43m, 
Wayne Rasband National Institute of Health, USA) in 2010-LD 
and 2010-HD). For 2010-LD and 2010-HD, we also determined 
the area of each leaflet by image analysis. For lateral leaflets, there 
was no significant difference between opposite leaflets for length, 
width or area, thus, for data analysis we used the mean value 
of length, width or area of each pair of opposite leaflets of an 
individual leaf. 

Time-course measurements of visible leaf number, internode and 
terminal leaflet lengths 

We randomly selected 33, 28, and 22 plants from the 2007-Su, 
2010-LD, and 2010-HD crops, respectively. These plants were not 
those used for destructive measurements. The time of bud break 
of the cutting (referred to simply as bud break hereafter), mark- 
ing the start of primary shoot development, was noted for each 
plant. Five times per week for 2007-Su and six times per week for 
2010-LD and 2010-HD, we counted the number of visible leaves 
on each plant, including the scaly and stipular leaves located at 
the base of the shoot. A leaf was counted as visible as soon as its 
tip emerged. On primary shoots, we measured the length of the 
terminal leaflet of each leaf and the length of each internode with 
a ruler, from the first day on which the organ was fully visible 
until the end of organ expansion. A terminal leaflet was fully vis- 
ible when its insertion on the rachis was visible, and an internode 
was fully visible when the node at its base was visible. Phytomer 
position of the organ was specified. The total number N p of phy- 
tomers, including the peduncle, of the primary shoot of the plant 
p was determined at the end of shoot expansion. 

DATA ANALYSIS 

Convention for phytomer position 

For data analysis, organ position along the shoot was expressed in 
terms of the relative rank of the phytomer (i), such that the most 
basal or apical phytomers could be compared regardless of the 
number of phytomers per shoot. Relative rank was calculated as: 

i=(r-l)/(N p -l) (1) 



Days after 
bud break 

15 72.3 ±3.1 8.5 ±0.2 

30 93.9 ± 1.0 37.8 ±2.6 

36 95.0 ±0.1 47.0 ±2.1 



20.2 ±0.9 
22.3 ± 1.6 22.8 ±1.4 



56.1 ±7.7 



where r is the absolute rank numbered from the base (r = 1) to 
the top (r = Np) of the shoot. The relative rank of the peduncle is 
therefore 1. For the figures, the rounded phytomer relative rank, 
the relative rank of the phytomer rounded to the nearest tenth 
(0.1), was used for the sake of clarity. A phytomer consists of an 
internode, the leaf at the top of the internode and its axillary bud. 

Time of leaf appearance 

We calculated thermal time from air temperature above the crop, 
using a base temperature of 2.1°C. This base temperature was 
determined for Rosa hybrida "Radrazz" in a previous experiment 
(Guerin, pers. com.) according to Yang et al. (1995). All variables 
defining the time of an event are expressed in degree days (°Cd) 
since cutting bud break. 

The thermal time of leaf appearance t p (i) for relative rank i of 
plant p could not be directly assessed by measurement, because 
the times at which the measurements were made did not cor- 
respond exactly to the time of leaf appearance. Thus, time t p (i) 
was estimated from the parameters of a piecewise linear func- 
tion fitted to the time of the first observation of each leaf (with 
or without leaflets) of plant p (Figure 1): 

[Up + fipip + y p (i - c p ) ifc p <!<l 

where c p is the relative rank at which the two lines intersect, 
defining the transition between a phase of rapid leaf appear- 
ance and a phase of slow appearance, a p is the thermal time 
at which the first leaf appeared, f$p and y p the phyllochrons 
(°Cd) for the leaves of relative ranks below and above c p , respec- 
tively. Peduncles do not bear leaves, but a date of "virtual 
leaf" appearance was calculated by linear extrapolation for the 
peduncle (z = 1). 

Fitting of organ expansion 

Each organ on the shoot was characterized by its kind (k, with k = 
lea for terminal leaflets and k = int for internodes) and relative 
rank i. If the final organ length exceeded 12 mm, the time course 
of length Lp_k(t, t), expressed in mm, was fitted with a logistic 
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0.4 0.6 0.8 

Relative phytomer rank 



FIGURE 1 | Time (in degree days since budbreak) at which a leaf was 
observed for the first time (symbols) plotted against relative phytomer 
rank for an individual plant. The lines correspond to the fitting of a 
piecewise linear function, as in Equation (2), to the experimental data. 



function of the thermal time (f ) since bud break of plant p: 



The reasons for using such polynomials are given in the results 
section. 

MODEL DESCRIPTION 

Overview 

The data analysis was used to develop a model reproducing the 
organ expansion kinetics of individual plants from a real crop of 
rose bush primary shoots, from input data. Organ length was sim- 
ulated as a logistic function of thermal time since bud beak. In 
the reference scenario (So), the input variables were chosen so as 
to capture variability in both final length and the time of organ 
expansion. Two other scenarios (Si and S2), differing in the num- 
ber of input variables, were compared with the reference scenario. 
In these three scenarios, the input variables selected were both 
easy to measure and non-destructive. 

Reference scenario So 

Scenario inputs. The input variables for each plant p are: 



L p , k(i, t) 



1 + expl 4 



a) 



(3) 



where L™ k (i) is the maximal length of the organ, t~ k (i) is thermal 
time at the inflexion point, referred to as time at mid-expansion 
hereafter, and is the expansion rate (mm °Cd _1 ) at the 

inflexion point (maximal expansion rate). The ratio w™j.(0 = 
v p" k^l jfc(0 i s m e maximal expansion rate when organ size 
is normalized, i.e., the maximal expansion rate of the ratio 
Lp, /t(i, O/^j-O)- ^ nna l organ length was lower than 12 mm, 
which includes the internodes at the very base of the shoot for 
most plants and the internode just below the peduncle for a few 
plants, the thermal time course of organ length was not fitted with 
any function, because of a lack of accuracy, and the organ was 
discarded from the analyses using the fitting data. 

For each individual organ whose final length exceeded 12 mm, 
we used the values of the parameters of the logistic function to cal- 
culate the duration of expansion, defined as the time required for 
the organ to expand from 10 to 90% of its final length. Duration 
of expansion of an organ is inversely proportional to the value of 
w™ k (i) for this organ. 



(11) The number N p of phytomers of p from which we 
derived the set of the relative ranks 0, l/(N p — 1), ... , 
{N p - 2)/(N p - 1), 1, using Equation (1). 

(12) The time tp'° hs (i) when the leaf of relative rank i (i = 
0, . . . , (Np — 2)/(Np — 1)) ofp was observed for the first time. 

(13) The final length L™'y s (i) of the terminal leaflet of relative 



rank i (i = 0, . . . , (Np 
(14) The final length L 
i(i = 0, . ..,l)ofp 



2)/(N p - 1)) ofp. 



. obi 

p, int 



(z) of the internode of relative rank 



Scenario outputs. The outputs of the model consist of predicted 
functions. For each plant p of Np phytomers, the predicted length 
of its organ k (k = lea for terminal leaflet and k = int for intern- 
ode) at relative rank i, expressed in thermal time t since bud break 
ofp is 



T m, obs / -n 

V* W 



l + exp(4<(0 (t° p r\r)-t)) 



(5) 



where the parameters f° '^ red (z') and w k (i) are specified in the next 
paragraph. 



Modeling time of the inflexion point of the expansion function 

The time f° k (i) of the inflexion point of the time course func- 
tion of the organ k (lea for terminal leaflet, int for internode) 
at relative rank z of plant p was estimated from observed data 
as described in Fitting of Organ Expansion. We also aimed to 
predict the time of the inflexion point c k (i) from the time of 
leaf appearance t£(i), therefore k (i) was calculated as polyno- 
mials of the time of leaf appearance ff(i) of the corresponding 
phytomer: 



f p, lea® = Sfea . 0 + Sfe ". 1 + 2 



,0 



Sim. 0 + $int, lfpO) 



. if i < 1 

fp ' mt I&mo + Smi^© if»=i 



(4.1) 
(4.2) 



Parameter computations. From the observed times fp' obs (z), we 
estimated the parameters dp, fi p , y p , and c p of Equation (2), 

then we deduced the adjusted times (tp' a of leaf appearance. 
These appearance times were used to compute the predicted times 
tp' V k ed (i) corresponding to the inflexion points, using Equations 
(4.1) and (4.2). 

First, we estimated the parameters w™ k (i) 

fitting procedure introduced in the paragraph Fitting of Organ 
Expansion for each organ kind, relative rank, and plant. Then, 
using these estimations for all plants, we estimated w™ (z) for each 
relative rank z with the LOESS method of SAS. In the model, the 
estimated value of w k '(i) thus depends on kind of organ and its 
relative rank, but not on the individual plant. 



tS-ttt using the 

p.k w 
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Scenario 1 

Scenario input. The time of leaf appearance is the most demand- 
ing of the input variables to measure, because it requires a series 
of observations to be taken at particular times, with little flexi- 
bility in timing possible. We therefore, decided to eliminate this 
observation from scenario Si . So the input variables of each plant 
p are (II), (13), and (14). 

Scenario outputs. The output functions of the scenario Si 

^m,obs 

^pred .. v p,k ^ J , . 

>• * U ' t] ~ l+exp(4<(0(l»(i)-t)) W 

are similar to the outputs of the scenario So Equation (5), except 
that the thermal time at mid-expansion is independent of the 
plant and depends only on relative rank. The parameter w™(z) 
is specified in the paragraph Reference Scenario So Parameter 
Computations and the parameter t®(i) is specified in the next 
paragraph. 

Parameter computations. In Si, the value fj?(i) was estimated by 
relating the values f° k (i), estimated with the fitting procedure 
introduced in the paragraph Fitting of Organ Expansion, with 
relative rank i by the LOESS method of SAS. 

Scenario 2 

Scenario input. Scenario S2 simulates the kinetics of expansion of 
plants differing only in phytomer number. The comparison of So 
with S2 was designed to estimate the gain in accuracy with the use 
of a model using n different individual plants to simulate expan- 
sion (So) rather than a model reproducing n mean plants differing 
only in phytomer number (S2). In S2, the only input variable is 
the number of phytomers per primary shoot (II), which is used 
to calculate relative ranks. 

Scenario output. The output of the scenario S2 is independent of 
the plant in the sense that the organ lengths of plants with the 
same number of phytomers are similar: 

k ( ' ) l+exp(4wr(0( t °(i)- t )) 

where the parameters w™(z) and f?(z) are previously specified and 
the parameter L™(0 is specified in the next paragraph. 

Parameter computations. The final organ length L™(i) was esti- 
mated by relating the measured values of L™ k (i) to i, with the 
LOESS method of SAS. 

MODEL EVALUATION 

The three scenarios were evaluated by the leave-one-out cross- 
validation method (Linhart and Zucchini, 1986). 

For each plant p of the original data set consisting of n plants: 

- We removed plant p and estimated the model parameters from 
the n — 1 remaining plants. 



- We computed the expansion curves for each terminal leaflet 
and each internode of plant p, using the scenario output, and 
calculated the mean square error of prediction: 

E?=f (^ k (i,t)-L^(i,t)) 2 

MSEP p . k (i) = ^ — -± '- (8) 

np,k(i) 

where n Pz ^(i) is the number of observations for an organ of kind k 
of plant p of relative rank i, i° b £(£, i) is the observed organ length 

at time t and (i, t) (or L^ red (i, t ) for S2 ) is the predicted organ 
length at time t computed from one of the functions given in 
Equations (5-7). For MSEPp k(i) calculations, we retained only 
four values for the plateau, defined as the time at which organ 
length exceeded 0.97 times the final length, to avoid giving too 
much weight to the plateau. 

Mean MSEPp^(i) was then calculated for all plants, both 
per rounded relative rank and for all phytomer relative ranks 
pooled together. The square root of MSEP (RMSEP) was then 
calculated. 

STATISTICAL ANALYSES 

Principal component analyses were performed with SPAD soft- 
ware (V7.4, Coheris, Suresnes, France). All other statistical anal- 
yses were performed with SAS software (V 9.3, SAS Institute). 
For the fitting of logistic functions Equation (3) and of piece- 
wise linear functions Equation (2), we used the NLIN procedure. 
For the establishment of relationships between relative rank i on 
the one hand and w™(i), t®(i), or L^ii) on the other, we used 
non-parametric methods and the LOESS procedure. 

RESULTS 

We will first present the variability present in the datasets used 
for expansion studies. We will then propose a simplified rep- 
resentation of rose compound leaves for studies of the kinet- 
ics of expansion, and describe the principal features of organ 
expansion kinetics: global pattern, origin and amplitude of the 
variability and how this variability can be captured through 
relationships. Finally, we will present a model based on this anal- 
ysis that simulates leaflet and internode expansion kinetics from 
non-destructive input data. 

ARCHITECTURAL VARIABLES WERE HIGHLY VARIABLE BOTH 
BETWEEN PHYTOMER POSITIONS AND BETWEEN PLANTS 

High levels of variation were observed in all three crops used for 
time-course measurements. Primary shoots comprised between 
10 and 16 phytomers on the plants used for time-course measure- 
ments (Table 2). The rate of leaf appearance was highly variable, 
with the phyllochron during the phase of slow leaf appearance 
varying from 20.6 to 53.4° Cd between plants within a crop, result- 
ing in considerable variation of the dates of leaf appearance at 
equivalent phytomer positions (e.g., at rank 10, leaf appearance 
on the first plant occurred 172°Cd before that on the last plant 
within a crop). Final lengths were between 6.5 and 87.5 mm for 
terminal leaflets and between 1.0 and 63.5 mm for internodes, 
considering all ranks together. At a given phytomer position, final 
lengths also varied considerably between plants. For example, 
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Table 2 | Mean values and variability of architectural variables for the three crops used for the expansion study. 





Phytomers 


Phyllochron during 


Time of leaf 


Final lengths (mm). 


Final lengths (mm) at 




per shoot 


the phase of slow 


appearance at rank 


all ranks together 


relative rank 0.8 






leaf appearance 


10 (°Cd since bud 














(°Cd) 


break) 


Terminal leaflet 


Internode 


Terminal leaflet 


Internode 






Mean ± std 


12.0 ± 0.9 c a 


35.0 ± 7.8 b 


198.0 ±28.3 a 


60.6 ± 12.2 ab 


21.0 ± 12.8 b 


71.5 ± 5.4 b 


27.8 ± 6.7 b 


Min-max 


10-13 


21.2-50.7 


152.3-264.3 


18.0-83.0 


1.0-62.0 


58.0-80.0 


9.0-40.0 




Mean ± std 


14.4 ± 1.0 a 


35.5 ± 5.9 b 


127.1 ± 24.4 c 


62.4 ± 16.7 a 


23.9 ± 14.1 a 


77.5 ± 5.1 a 


36.0 ± 5.5 a 


Min-max 


13-16 


24.5-46.8 


79.0-173.4 


6.5-87.5 


3.5-61.5 


63.0-86.0 


22.5-48.0 




Mean ± std 


13.5 ± 1.4 b 


39.7 ± 7.9 a 


159.5 ± 47.6 b 


58.2 ± 16.0 b 


24.5 ± 16.0 a 


71.0 ± 6.2 b 


34.9 ± 7.7 a 


Min-max 


10-16 


20.6-53.4 


87.0-258.6 


6.5-82.0 


3.0-63.5 


59.0-82.0 


16.0^9.5 



"Mean values with the same letters do not differ significantly between crops at P < 0.05, in One-Way ANOVA followed by LSD tests. Variability was assessed by 
determining both the standard deviation and the minimal and maximal values within each crop. 



at a rounded relative rank of 0.8, the range of final lengths 
was 58.0-86.0 mm for terminal leaflets and 9.0-49.5 mm for 
internodes. Differences in mean values between crops, for phy- 
tomer numbers, timing of leaf appearance and final lengths, were 
significant. 

TERMINAL LEAFLET LENGTH CAN BE USED TO SIMULATE THE 
DIMENSIONS OF ALL LEAFLETS WITHIN A LEAF 

Rose leaves have generally odd numbers of leaflets, with between 
one and nine leaflets per leaf, except at the base of the shoot, 
where leaves are reduced to scales or stipules. Principal compo- 
nent analysis (PCA) was performed on the lengths, widths and 
square roots of the areas of all individual leaflets of both grow- 
ing and mature leaves collected from the 2007-Sp, 2010-LD, and 
2010-HD crops. The percentage of the variance explained by the 
first principal component was 97, 79, 79, and 80% for 1-, 3-, 5-, 
and 7-leaflet leaves, respectively (no nine-leaflet leaf sampled in 
these experiments). These very high values for a single compo- 
nent indicate that all the dimensions (lengths, widths and square 
roots of areas) of all the leaflets were correlated within a given 
leaf, for both growing and mature leaves. Consequently, we were 
able to simplify the study of leaf expansion by focusing on a single 
dimension. We chose to study terminal leaflet length, because all 
leaves composed of leaflets have a terminal leaflet and this dimen- 
sion is easily measured earlier in expansion than the other possible 
dimensions. 

Terminal leaflet length (L) can be used to calculate whole leaf 
area (A) at any growth stage, given the number of leaflets of the 
leaf (Ni), from the following equation: 

A = 0.287 L 2 N° 746 (n = 201, R 2 = 0.93) (9) 

With a multiplicative model for terminal leaflet length and leaflet 
number per leaf, the residual distribution did not depend on 
leaflet number per leaf as with a linear model. The two parameter 
values (0.287 and 0.746) were adjusted so as to maximize R 2 . 
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FIGURE 2 | Thermal time course of terminal leaflet (diamond) or 


internode (triangle) length for an individual organ. Symbols represent 


experimental data and the lines correspond to the symmetric logistic 


function fitted to the experimental points. 





EXPANSION FOLLOWS A SIGM0IDAL PATTERN THAT CAN BE FITTED 
BY LOGISTIC FUNCTION 

In individual organs, changes in terminal leaflet or internode 
length as a function of thermal time followed a sigmoidal pat- 
tern (Figure 2). This pattern of change was well fitted by a 
logistic function (Table 3): R 2 close to 1 and RMSEP around 
1 mm. For all phytomers, the period of internode expansion was 
included within the expansion period of the terminal leaflet of 
the same phytomer (Figure 2). Terminal leaflets reached 10% 
of their final size a mean of 31°Cd before the correspond- 
ing internodes reached 10% of their final size, and the leaflets 
reached 90% of their final size a mean of 35°Cd after the 
internodes. 

DIFFERENCES IN EXPANSION KINETICS ARE LARGELY DUE TO 
DIFFERENCES IN THE TIME AT WHICH EXPANSION OCCURS AND TO 
MAXIMAL EXPANSION RATE 
Between plants at a given phytomer position 

Within each relative rank, the variability of terminal leaflet 
expansion kinetics was high, as shown for relative rank 0.7 
in Figure 3A. This variability resulted from high variability of 
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Table 3 | Goodness-of-fit of logistic functions to expansion curves of 
individual terminal leaflets and internodes: number of fitted curves, 
mean values of Ft 2 , root mean square error (RMSE), and mean 
number of observations per curve. 



Terminal leaflets Internodes 



Number of organs for which expansion 

curves were fitted 

ft 2 

RMSE (mm) 

Number of observations per curve 



622 

0.99 

1.2 

19 



711 

0.98 
0.81 
16 
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FIGURE 3 | Elongation of individual organs positioned at relative 
rounded rank 0.7, as a function of thermal time since bud break, for 
terminal leaflets (A) or internodes (C). (B,D) show the same data for 
normalized length (length at any time divided by final length) and for 
thermal time counted since the inflexion point of the expansion curve of 
each individual organ. The crops 2007-Su, 2010-LD, and 2010 HD are 
represented by different symbols and colors. Each point corresponds to the 
measurement of an individual organ. Examples of individual curves are 
given for three plants (plants 11, 34, and 76). 



both size and the time at which expansion occurred (265° Cd 
separating the times at mid- expansion of the first and last 
plants). All individual observations were almost superimpos- 
able if the length of each organ at a given time was nor- 
malized by dividing by its final length, and if the time was 
expressed as time since mid-expansion of individual organ 
(tp k (i); Figure 3B). Thus, the differences in organ size between 
plants essentially resulted from differences in maximal expansion 
rate, rather than differences in expansion duration. This was true 
for all relative ranks, for both terminal leaflets and internodes 
(Figures 3C,D). 

Between phytomer positions 

For terminal leaflets, the variations of final length and maxi- 
mal expansion rate with phytomer position followed the same 
tendencies for all three crops: a strong increase in the basal 
part of the shoot (up to relative ranks 0.4-0.5), then a more 
moderate increase up to relative rank 0.8, followed by a 
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FIGURE 4 | Final length (A,C), maximum expansion rate (B,D, closed 
symbols) and expansion duration (B,D, open symbols) plotted against 
rounded relative phytomer rank for terminal leaflets (A,B) and 
internodes (CD). The 2007-Su, 2010-LD, and 2010-HD crops are 
represented by different symbols and colors: black triangles for 2007-Su, 
blue circles for 2010-LD and red squares for 2010-HD. Error bars represent 
the standard deviation, with 33, 28, and 22 plants in the 2007-Su, 2010-LD, 
and 2010-HD crops, respectively. 



slight decrease or stability (Figures 4A,B). Expansion duration 
increased moderately from the base to relative rank 0.4 (2007-Su) 
or 0.7 (2010-LD and 2010-HD) and then stabilized and decreased 
for upper ranks (Figure 4B). 

For internodes, final length (Figure 4C) increased up to rel- 
ative ranks 0.8 then strongly decreased for relative rank 0.9 and 
increased strongly for relative rank 1.0 (peduncle). The variations 
of maximal expansion rate followed an essentially similar pat- 
tern to final length, except at the base of the axis (Figure 4D). 
By contrast, variations of expansion duration along the shoot 
(Figure 4D) were characterized by stability from relative ranks 
0.3-0.9, with a lower duration at relative rank 0.2 and a much 
greater duration (x 1.7) for the peduncle. 

For both terminal leaflets and internodes, absolute values of 
final length, maximal expansion rate and expansion duration dif- 
fered between the crops; in addition, the ranking of the three 
crops for these variables differed between different zones of the 
shoot. However, for all crops, the amplitude of variation between 
phytomer relative ranks was much higher for final length and 
maximal expansion rate than for expansion duration. 

Within each crop, plants with a high leaflet or internode final 
length at a particular phytomer position did not necessarily also 
have a greater length at the next phytomer (not shown); the same 
was true for low leaflet or internode final lengths. 

TIME OF LEAF APPEARANCE CAN BE USED TO ESTIMATE TIME AT 
MID-EXPANSION 

Figure 5 shows the mid-expansion times for individual leaflets 
(n = 629) and internodes (n = 584) plotted against time of leaf 
appearance. These relationships integrate both the rate at which 
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expansion occurs on successive phytomers and the duration of 
expansion, because the predicted variable is the time at the middle 
and not at the beginning of expansion. For terminal leaflets, 
the relationship between time at mid-expansion and time of 
leaf appearance gradually curves, because expansion duration 
changes gradually with phytomer relative rank (Figure 4B). We 
therefore decided to use second-order polynomials to adjust the 
relationship for terminal leaflets. For the internodes, expansion 
duration was almost stable for vegetative internodes and much 
longer for peduncles (Figure 4D). We therefore, used first-order 
polynomials to adjust the relationship between time at mid- 
expansion and time of leaf appearance, using different parameters 
for vegetative internodes and peduncles. Estimates for the param- 
eters of the relationships defined in the Equations (4.1) and (4.2) 
are given in Table 4. 

The intercepts 8;^ 0 and &,„ f o correspond, for phytomers 
whose leaf appears at bud break, to the time lag between leaf 
appearance and organ (leaflet or internode, respectively) maxi- 
mal expansion rate. The other parameters of Equations (4.1) and 
(4.2) have no direct biological meaning. The adjusted determina- 

0.97 (n = 627) for terminal leaflets, 

525) for vegetative internodes and R^. = 0.86 
(n = 59) for peduncles. These high R 2 indicate that time of leaf 
appearance can be used to estimate time at mid-expansion for 
individual leaflets and internodes, because these relationships 



tion coefficients were B~^. 



capture the variability between crops, between individual plants 
and between phytomer positions. 

A MODEL OF LEAFLET AND INTERNODE EXPANSION REPRODUCES 
INTERPLANT VARIABILITY 

From the results presented above, we developed a model repro- 
ducing the terminal leaflet and internode expansion of a crop 
of rose bush primary shoots. The three scenarios of the model, 
differing in the number of input variables, were evaluated by 
cross-validation. 

In the reference scenario (So), input variables accounted for 
interplant variability in both organ final lengths and time of leaf 
appearance. The term w™(i) depends only on phytomer posi- 
tion and does not vary between individual plants. This term 
represents the maximal expansion rate for a normalized expan- 
sion curve and is inversely proportional to expansion duration. 
The parameters of So are given in Table 4 and the values of 
w™(i) estimated by the LOESS method are given in Table 5. So 
correctly reproduced the expansion kinetics of the different termi- 
nal leaflets and internodes of an individual plant (Figures 6A,C) 
and of different plants within a given crop (Figures 6B,D). The 
RMSEP for all ranks together was 4.5 mm for leaflets and 2.5 mm 
for internodes, corresponding to 7.3 and 10.2%, respectively, of 
organ final lengths. The accuracy of the model was similar for the 
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FIGURE 5 | Time at the inflexion point of the expansion curve for 
leaflets (A) or internodes (B), plotted against the time of leaf 
appearance. All times are expressed in degree days since bud break. Each 
point corresponds to an individual organ. Data corresponding to the same 
rounded relative phytomer are represented by the same symbol, without 
distinction between the three crops (2007-Su, 2010-LD, and 2010-HD). The 
lines correspond to Equation (4), parameter values are given in Table 5. 



Table 4 | Values of the model parameters used in the reference 
scenario (So) other than those obtained by the LOESS method (see 
Table 5). 



Parameter 

8/e3,0 (°Cd) 

5 fea,1 

5,„ a , ("Cd" 1 ) 



:°cd) 



°fea,2 
8 /nt,0 
$;nf,1 

8ped,0 (°Cd) 
5ped,1 



Equation 

4.1 
4.1 
4.1 
4.2 
4.2 
4.3 
4.3 



Value 

61.5 
1.62 
-0.00113 
70.6 
1.34 
220 
0.893 



Table 5 | Model parameters: values of w™(i), t°(i), and L?(i) 
estimated from the experimental data by the LOESS method. 



Relative 



Terminal leaflets 



Internodes 



ranK 


W,™ (/) X 10 3 

<°Cd) 


(°Cd) 


(mm) 


WP nt (i) x 10 3 
(°Cd) 


Co 

(°Cd) 


C t C) 

(mm) 


0.0 












5.9 


0.05 












7.0 


0.1 










23.7 


9.0 


0.15 


7.39 


74.7 


19.6 




38.8 


10.4 


0.2 


7.16 


82.4 


27.1 


11.72 


54.9 


11.7 


0.25 


6.94 


89.7 


34.3 


10.99 


70.5 


12.5 


0.3 


6.72 


97.1 


41.5 


10.35 


86.1 


13.0 


0.35 


6.53 


105.0 


46.9 


9.76 


100.4 


13.3 


0.4 


6.31 


116.6 


52.5 


9.27 


118.8 


14.6 


0.45 


6.11 


132.9 


55.3 


8.88 


138.8 


16.1 


0.5 


5.93 


154.3 


57.8 


8.66 


162.6 


17.9 


0.55 


5.78 


180.7 


60.1 


8.76 


189.4 


20.4 


0.6 


5.66 


211.6 


62.7 


8.86 


217.7 


23.4 


0.65 


5.56 


250.3 


65.9 


8.95 


250.3 


28.2 


0.7 


5.53 


281.9 


68.4 


8.97 


277.0 


31.7 


0.75 


5.54 


313.6 


71.1 


8.94 


306.1 


33.6 


0.8 


5.60 


340.2 


73.1 


8.77 


334.8 


32.4 


0.85 


5.70 


362.6 


73.4 


8.43 


363.9 


26.1 


0.9 


5.79 


381.4 


72.2 


7.94 


394.6 


14.4 


0.93 


5.84 


393.3 


71.4 


7.34 


414.2 


8.0 


1.0 








5.19 


457.6 


52.7 



w^(i) is used in scenarios So, S; and S2; ti(i) is used in scenarios S7 and S?; 
UFO) is used in scenario S?. 



Frontiers in Plant Science | Plant Biophysics and Modeling 



October 2013 | Volume 4 | Article 418 | 8 



Demotes-Mainard et al. 



Organ expansion in rose bush 



different phytomer positions, for both leaflets (Figure 7A) and 
internodes (Figure 7B). 

Scenario 1 was designed to assess model accuracy assuming 
that the time at which expansion occurred depended only on 
phytomer position, rather than being estimated, as in So, from 
an input variable. The parameters of Si are given in Table 5. 
In Si, the RMSEP for all ranks together was 8.6mm for leaflets 
and 4.5 mm for internodes (14.2 and 18.2% of final lengths, 
respectively). The RMSEP calculated for each phytomer position 
increased from the base to the top of the shoot (Figure 7). This 
was due to an increase in variability between plants for mid- 
expansion time, from the base to the top of the shoot (Figure 5). 
The optimal scenario would thus involve estimating the time at 
mid-expansion from a value dependent only on relative rank for 
basal phytomer positions and from measurements of leaf appear- 
ance time for more central and apical positions. In practical 
terms, this would decrease measurement time with respect to So, 
because the leaf-like structures at the base of the shoot are the 
most difficult to observe. 



Scenario 2 simulates the kinetics of expansion of plants dif- 
fering only in terms of phytomer number. The parameters of S2 
are given in Table 5. For S2, RMSEP was calculated by plant, so 
the comparison of the RMSEP values for So and S2 provides an 
estimate of the gain in accuracy when a scenario accounting for 
interplant variability (So) is used to simulate the expansion kinet- 
ics of a heterogeneous crop over a scenario reproducing n mean 
plants. In S 2 , RMSEP for all ranks together was 10.6 mm (17.4 %) 
for leaflets and 6.0 mm (24.2%) for internodes, these values being 
2.4 times higher than those for So- This loss of accuracy for S2 
with respect to So concerned all phytomer positions (Figure 7). 
It should be borne in mind that the method used to calculate 
MSEP does not aim to assess the ability of S2 to simulate a mean 
plant, because simulated lengths were not compared with mean 
observed lengths. 

DISCUSSION 

The aim of this study was to assess the kinetics of expansion of 
rose bush primary shoot internodes and leaves and to propose an 
empirical model reproducing organ expansion kinetics from non- 
destructive, easy-to-measure variables. Bush roses are ornamental 
and are sold as individual plants. We therefore, paid particular 
attention to interplant variability in expansion kinetics and to the 
ability of the model to simulate this variability. 

GROWING CONDITIONS AND INTER- AND INTRA-CROP VARIABILITY 

We studied expansion in three crops (2007-Su, 2010-LD, and 
2010-HD) subjected to different growing conditions, contrasting 
in terms of incident PAR, humidity and plant density in particu- 
lar, to ensure that we obtained robust results. Contrasted exper- 
imental conditions were used to find regularities in expansion 
kinetics, but not to establish quantitative relationships between 
architectural traits and environmental factors, such as density, 
light intensity or quality. This would require specific experiments. 
The large differences in PAR interception between densities in 
2010 (Table 1) suggest that light quality in the canopies, including 
the red-far red ratio in particular, differs between plant densities. 
These differences in growing conditions resulted in significant 
differences between crops in terms of the mean values of archi- 
tectural variables. Within each growing condition, the plants also 
displayed considerable variability. This level of interplant variabil- 
ity is typical of rose bush crops. This variability was observed 
although we selected cuttings to reduce heterogeneity due to 
topophysis (Bredmose et al, 2001) and differences in cutting leaf 
area (Costa and Heuvelink, 2003). 

AT ANY TIME DURING EXPANSION, THE DIMENSIONS OF THE WHOLE 
LEAF CAN BE RECONSTITUTED FROM TERMINAL LEAFLET LENGTH 

Our results show that it is possible to reconstitute the dimensions 
of the different leaflets of a rose leaf from a single dimension at any 
time during leaf expansion, given the number of leaflets per leaf. 
These results generalize those of Gao et al. (2012), who established 
that, for fully expanded leaves, total leaf area could be inferred 
from one leaf dimension and number of leaflets, within a geno- 
type. The allometric relationships using terminal leaflet length 
as a predictor of either whole leaf area (presented here, section 
Terminal Leaflet Length can be Used to Simulate the Dimensions 
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FIGURE 6 I Observed (symbols) and simulated (lines) organ lengths as 
a function of thermal time for the different relative ranks of an 
individual plant (A,C) or for different plants at the same relative rank 
(B,D). Illustrations (A,B) show terminal leaflet lengths, (CD) show 
internode lengths. As an example, an individual plant was chosen for 
illustrations (A,C), and a rounded relative rank of 0.7 is used for illustrations 
(B,D). Simulated values were obtained with the cross-validation method. 
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FIGURE 7 I Square root of the mean square error of prediction 
(RMSEP) for terminal leaflet (A) and internode (B) expansion, plotted 
against rounded relative phytomer rank. The short horizontal lines show 
the value of RMSEP for all phytomer positions considered together. RMSEP 
was calculated throughout the expansion of each individual organ, 
according to the three scenarios of the model. In So (reference scenario) 
input variables account for the interplant variability in both organ final length 
and time of leaf appearance; in Si , input variables account for the interplant 
variability in organ final length; S2 simulates mean plants. 
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of all Leaflets Within a Leaf) or leaflet dimensions (Demotes- 
Mainard et al., 2009) should probably be established for each 
genotype. Indeed, only relationships based on the product of 
length by width have been found to be stable across a range of 
genotypes, both at the leaflet level (Rouphael et al., 2010) and the 
leaf level (Gao et al., 2012), because the shape of leaflets and their 
insertion position along the rachis (loose or close) vary between 
genotypes. 

STABILITY OF ORGAN EXPANSION DURATION AND VARIABILITY OF 
FINAL ORGAN SIZE 

Differences in organ final size between phytomer positions and 
between plants at a given position resulted mainly from differ- 
ences in maximal expansion rate, rather than expansion duration. 
This is similar to the behavior of other species. For example, 
whereas final leaf and internode size widely varies between phy- 
tomers along the main shoot, organ expansion duration little 
varies in narrow-leafed lupin (Dracup and Kirby, 1993) and is 
stable in arabidopsis (Mtindermann et al., 2005). At a given phy- 
tomer position, expansion duration is stable in thermal time 
between plants grown in contrasted growing conditions for 
sunflower leaves (Granier and Tardieu, 1998; Dosio et al., 2003), 
sorghum leaves (Lafarge and Tardieu, 2002) and internodes (Xue 
et al, 2012). In tall fescue and wheat, expansion duration is more 
stable if expressed in phyllochronic time (thermal time divided 
by the phyllochron) than if expressed in thermal time (Fournier 
et al., 2005). This stability in phyllochronic time has been inter- 
preted as an emerging property of a self-regulated system, in 
which the appearance of an organ triggers changes in expansion 
(Fournier et al., 2005; Verdenal et al, 2008). Under this hypoth- 
esis, the link between phyllochron and duration results from a 
large proportion of expansion occurring when the leaf is still in 
the whorl generated by previous leaves. 

In rose bush, as in many other species (for example lupin, 
Dracup and Kirby, 1993; wheat, Evers et al., 2005; cucumber, 
Kahlen, 2006; sorghum, Xue et al, 2012), there was a general gra- 
dient of final organ length with phytomer relative rank common 
to all three crops and plants. This gradient can be used to model 
a mean plant, as in scenario S2. However, this gradient varied 
between individual plants: a plant with a high (or low) leaflet 
or internode final length at a particular phytomer position did 
not systematically present a high (or low) length at the next phy- 
tomer. This was unexpected, at least for basal phytomer positions. 
Indeed, during their expansion basal phytomers depend on the 
cutting, both for nitrogen, which is provided by remobilization 
(Cabrera, 2003) and for photosynthesis (Costa and Heuvelink, 
2003). In addition, on average eight phytomers closest to the 
base of the plant are preformed in the bud in the variety studied 
here (Girault et al, 2008), suggesting that they may be influ- 
enced similarly by the physiological state of the mother stem. We 
therefore, assume that the differences in leaflet or internode final 
size between plants may reflect differences in the phylloclimate 
perceived by the organs during their expansion. Indeed, final size 
is dependent, in several species, on phylloclimate (or climate) 
at the time of organ expansion (Granier and Tardieu, 1999, for 
sunflower leaves; Gautier et al, 2000, for white clover internodes, 
petioles and leaves; Andrieu et al., 2004, for maize leaves; Kahlen 



and Sttitzel, 201 1, for cucumber internodes). Our model could be 
used to test this hypothesis, by accurately defining the timing of 
individual organ expansion. 

MODELING THE EXPANSION OF INDIVIDUAL PRIMARY SHOOTS 

We developed a model reproducing the kinetics of leaflet and 
internode expansion on rose bush primary shoots, based on non- 
destructive and easy-to-measure input variables. The comparison 
of the three scenarios showed that, to reproduce accurately the 
variability of expansion within a crop, it was necessary to mea- 
sure the final lengths of terminal leaflets and internodes and to 
count the number of phytomers per shoot and of leaflets per leaf, 
all of which can be done after the final primary shoot has ceased 
to elongate. The timing of leaf appearance must also be estab- 
lished, at least for leaves in the middle and upper part of the 
shoot, as fixed values of time at mid-expansion can be used for the 
basal phytomers. These data are easy to acquire, but their acquisi- 
tion is labor-intensive. The model based on these input variables 
satisfactorily simulated the diversity of individual plants. For eco- 
physiological studies, this model will make it possible to relate 
traits of interest measured on specific plants to features of expan- 
sion for the same plants, which is important for experimental 
crops with high levels of interplant variability. 

In order to represent architectural development, our model 
of organ expansion can be easily integrated in a FSPM of rose 
bush using L-systems. The addition of a three-dimensional struc- 
ture to the expansion model implies to model organs' shape 
and spatial position, as done for FSPMs of other species (e.g., 
Evers et al, 2005; Mtindermann et al., 2005; Kahlen et al., 
2008). This 3D model would provide essential information to 
study the environmental regulation of bud break, which strongly 
influences plant architecture and still remains a major research 
area (Domagalska and Leyser, 2011). Indeed, when coupled 
with a light model, the model will give information about 
light phylloclimate, which is a key factor controlling bud break 
through carbohydrate availability (Girault et al, 2010; Henry 
et al, 2011; Rabot et al., 2012), and light quality (Mor and 
Halevy, 1984; Girault et al., 2008). In rose bush, there is con- 
siderable interplant variability in the number and position of 
buds breaking along the primary shoot (for instance, in 2007- 
Su between 2 and 8 buds outgrew on primary shoots with the 
same number of axillary buds). It may therefore, be particu- 
larly relevant to be able to investigate the relationships between 
light phylloclimate and bud break at the level of individual 
plants. 

The differences in organ expansion kinetics between plants 
within a relative rank resulted mainly from differences in max- 
imal expansion rate, rather than differences in expansion dura- 
tion. This had two related implications, which were used in the 
model. Firstly, variations in organ final size were sufficient to 
capture interplant variability in size at any time during expan- 
sion, within a relative rank, as in scenarios So and Si. Secondly, 
in the model, the term w™(i), which is inversely proportional 
to expansion duration, can be predicted from a single value 
common to all plants within a relative rank with little loss of 
accuracy (all scenarios). In the model the term w™(i) depends 
on relative rank. However, the variations of expansion duration 
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between phytomer positions were moderate for terminal leaflets 
and for vegetative internodes. In scenario So, if w™(i) was 
replaced by only three values, one common to all terminal leaflets 
("Wy = 5.96 x 10~ 3 °Cd), one common to all vegetative intern- 



odes (w£ ( = 9.04 x 10 



-3°r 



Cd) and one for peduncles (w" ed = 

5.17 x 10~ 3 °Cd), the accuracy of the model remained the same 
(RMSEP = 7.4 and 10.2% of final length, for terminal leaflets and 
internodes, respectively, data not shown) as with w™(z). Thus, if 
the model has to be calibrated for a new genotype, this process 
could be simplified by estimating the model parameters from a 
limited number of phytomer positions. 

The time of leaf appearance has been used in several archi- 
tectural models for the prediction of ontogenic development (for 
example in wheat: Fournier et al., 2003; Evers et al., 2005; in 
cucumber: Kahlen, 2006; in ryegrass: Verdenal et al, 2008). In our 
model, we used this parameter in the reference scenario So to pre- 
dict the mid-expansion time for both leaflets and internodes. This 
coordinates expansion between different phytomer positions for 
the same kind of organ, and between leaflets and internodes. 

CONCLUSION 

This work provides insight into the kinetics of expansion at the 
level of individual organs that was lacking for rose bush. It is 
original in that it considers interplant variability, which has been 
little studied, in any plant species. This variability is important 



for rose bushes, because crops display high levels of interplant 
variability for architecture and the commercial product is the 
individual plant. On the basis of these results, we propose an 
empirical model that accurately reproduces, for individual plants, 
the expansion kinetics of primary shoots from non-destructive 
input variables. Even in its current state, this model already pro- 
vides a useful tool for studying ecophysiological processes, such 
as bud break response to light phylloclimate, because a primary 
shoot constitutes an interesting model of the whole plant, with 
simple interactions between organs. With a view to studying rose 
bush architecture, these results can be used as a grid for analyz- 
ing the expansion of shoots resulting from branching, and for 
investigating genotype differences. 
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